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Magnetic skyrmions are particle-like topological excitations in ferromagnets, which have the topologi¬ 
cal number Q = ± 1, and hence show the skyrmion Hall effect (SkHE) due to the Magnus force effect 
originating from the topology. Here, we propose the counterpart of the magnetic skyrmion in the anti¬ 
ferromagnetic (AFM) system, that is, the AFM skyrmion, which is topologically protected but without 
showing the SkHE. Two approaches for creating the AFM skyrmion have been described based on micro- 
magnetic lattice simulations: (i) by injecting a vertical spin-polarized current to a nanodisk with the AFM 
ground state; (ii) by converting an AFM domain-wall pair in a nano wire junction. It is demonstrated that 
the AFM skyrmion, driven by the spin-polarized current, can move straightly over long distance, benefit¬ 
ing from the absence of the SkHE. Our results will open a new strategy on designing the novel spintronic 
devices based on AFM materials. 


PACS numbers: 75.50.Ee, 75.78.Cd, 75.78.Fg, 12.39.Dc 


Skyrmion is a topologically protected soliton in continuous 
field theory, which is recently realized in both bulk non- 
centrosymmetric magnetic materials [1,2] and thin films [3], 
where the ferromagnetic (FM) background is described by the 
non-linear sigma model with the Dzyaloshinskii-Moriya in¬ 
teraction (DMI) [4]. The study of the magnetic skyrmion is 
one of the hottest topics in condensed matter physics, due to 
its potential applications in information processing and com¬ 
puting [5, 6]. There are several ways to create magnetic 
skyrmions, e.g., by applying spin-polarized current to a nan¬ 
odisk [7, 8], by applying the laser [9], from a notch [10] and 
by the conversion from a domain wall (DW) pair [11, 12]. A 
magnetic skyrmion can be driven by the spin-polarized cur¬ 
rent [13, 14]. However, it does not move parallel to the in¬ 
jected current due to the skyrmion Hall effect (SkHE), since 
its topological number is ±1. This will pose a severe chal¬ 
lenge for realistic applications which require a straight mo¬ 
tion of magnetic skyrmions along the direction of the applied 
current [13]. 

In this work, we demonstrate that a skyrmion can be nu¬ 
cleated in antiferromagnets, as illustrated in Fig. 1, based on 
micromagnetic lattice simulations. We refer to it as an antifer¬ 
romagnetic (AFM) skyrmion. We further show that the AFM 
skyrmion can move parallel to the applied current since the 
SkHE is completely suppressed, which is very promising for 
spintronic applications. 

Recently, antiferromagnets emerge as a new field of spin- 
tronics [15-18]. A one-dimensional topological soliton in 
antiferromagnets is an AFM DW [19]. An AFM DW can 
be moved by spin transfer torque (STT) induced by spin- 
polarized currents or spin waves [20, 21]. Besides, a two- 
dimensional (2D) topological soliton, that is, the magnetic 
vortex, has been studied in 2D AFM materials [22]. The AFM 
system has an intrinsic two-sublattice structure. The spins of 
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Figure 1. Illustration of an AFM skyrmion spin texture in a 
2D AFM monolayer with the two square antiparallel magnetic 
sublattices. The color scale represents the magnetization direction, 
which has been used throughout this paper: orange is into the plane, 
green is out of the plane, white is in-plane. 


the ground state are perfectly polarized in each sublattice. We 
may calculate the topological number of the spin texture pro¬ 
jected to each sublattice. Hence, we propose to assign a set of 
topological numbers (+1,-1) to one AFM skyrmion, which 
shows no SkHE since it has no net magnetization. 

We also present two approaches to create an AFM 
skyrmion. One is applying a spin-polarized current perpen¬ 
dicularly to a disk region, which flips the spin in the applied 
region. The other is a conversion from an AFM DW pair 
in junction geometry as in the case of the conversion of a 
FM skyrmion from a FM DW pair [11]. Furthermore, we 
show that it is possible to move an AFM skyrmion by apply¬ 
ing a spin-polarized current. The AFM skyrmion can travel 
very long distance without touching the sample edges. It is 
also insensible to the external magnetic field. These results 
will be important from the applied perspective of magnetic 
skyrmions. 
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Results 

AFM system. We investigate the AFM system with the lattice 
Hamiltonian, 


-Hafm = J y + (ttij xmj)-K ^(mf ) 2 , 

(ij) (i-J) * 

( 1 ) 

where rrii represents the local magnetic moment orientation 
normalized as \rrii\ = 1 , and (i,j) runs over all the nearest 
neighbor sites. The first term represents the AFM exchange 
interaction with the AFM exchange stiffness J > 0. The sec¬ 
ond term represents the DMI with the DMI vector D. The 
third term represents the perpendicular magnetic anisotropy 
(PMA) with the anisotropic constant K. 

The dynamics of the magnetization rrii is controlled by 
applying a spin current in the current-perpendicular-to-plane 
(CPP) configuration [14, 23]. We numerically solve the 
Landau-Lifshitz-Gilbert-Slonczewski (LLGS) equation, 


drrii 

dt 


= — \^\rrii x Hf f + oirrii x 


drrii 

dt 


+ | 7 | /3(m; xpx rrii), 


(2) 


where Hf r = — dH AFM /dm is the effective magnetic field 
induced by the Hamiltonian equation (1), 7 is the gyromag- 
netic ratio, a is the Gilbert damping coefficient originating 
from spin relaxation, (3 is the Slonczewski-like STT coef¬ 
ficient, and p represents the electron polarization direction. 
Here, (3 = | | with /i 0 the vacuum magnetic permit¬ 

tivity, d the film thickness, Ms the saturation magnetization, 
and j the current density. We take the p = — z for creat¬ 
ing the AFM skyrmion, while p = — y for moving the AFM 
skyrmion. Although an antiferromagnet comprises complex 
two sublattices of reversely-aligned spins, the STT can be 
applicable also for the AFM system provided the lattice dis¬ 
creteness effect is taken into account with an ultra-small mesh 
size in the micromagnetic simulations [17, 18]. The STT is 
induced either through spin-polarized current injection from 
a magnetic tunnel junction polarizer or by the spin Hall ef¬ 
fect [14, 24]. We can safely apply this equation for the AFM 
system since there is no spatial derivative terms. 

A comment is in order. We cannot straightforwardly use 
the current-in-plane (CIP) configuration to control the dynam¬ 
ics of the magnetization as it stands, since spatial derivative 
terms are involved in the LLGS equation, that is, urrii x 
(rrii x d x rrii ) + u'rrii x d x rrii. 


Topological stability. The skyrmion carries the topological 
number. In the continuum theory it is given by 

Q = [ d 2 xm(x ) • ( d x m(x ) x d y m(x )). (3) 

47 T J 

However, the AFM system has a two-sublattice structure 
made of the A and B sublattices. In our numerical compu¬ 
tation we employ the discretized version of the topological 
charge equation (3), 

Qr = m i ' ( m J X m fc) ’ (4) 

ijk 
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Figure 2. Creation of an isolated AFM skyrmion in a 2D AFM 
nanodisk using a vertically injected spin-polarized current, (a) 

AFM skyrmion in the nanodisk with a diameter of 80 nm created 
by a 2-ns-long spin-polarized current pulse (j = 5 x 10 13 A m -2 ) 
perpendicularly injected into the nanodisk in a circle region with dif¬ 
ferent size followed by a 1-ns-long relaxation (see Supplementary 
Movie 1). The spin-polarized current injection region is denoted by 
the blue circle, of which the diameter equals 30 nm, 40 nm, 50 nm, 
respectively. The size of all AFM skyrmions is found to be identical 
irrelevant of the current injection region, (b) Magnetization distri¬ 
bution of the AFM skyrmion in an AFM nanodisk. It is made of a 
toroidal DW with fixed radius and width determined by the material 
parameters. The AFM skyrmion size d s is defined by the diameter of 
the white circle, where m z = 0. (c) A 2-ns-long spin-polarized cur¬ 
rent pulse with different current density j is perpendicularly injected 
into the nanodisk with a diameter of 80 nm followed by a 1-ns-long 
relaxation. The initial state of the nanodisk is the AFM ground state. 
The spin-polarized current injection region is denoted by the blue 
circle, of which the diameter equals 40 nm. (d) Total micromagnetic 
energy E tota i for an isolated AFM skyrmion and the AFM ground 
state as a function of the DMI constant D. The D and E totSL \ are in 
unit of 10 -21 J and 10“ 17 J, respectively. 


for each sublattice (r = A, B). Hence, we propose to assign 
a set of two topological numbers (Qa, Qb) to one skyrmion. 
We obtain Qa = — Qb = 1 for a skyrmion in a sufficiently 
large area. Even if the skyrmion spin texture is deformed, its 
topological number cannot change. A skyrmion can be neither 
destroyed nor separated into pieces, that is, it is topologically 
protected. 

Creation of an AFM skyrmion by a vertical spin current. 

We employ a CPP injection with a circular geometry in a nan¬ 
odisk. The CPP injection induces spin flipping at the current- 
injected region. When we continue to apply the current, the 
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Figure 3. The time evolution of the topological number Q , that 
is, the skyrmion number, in the nucleation process of an isolated 
AFM skyrmion in a 2D AFM nanodisk with the two square sub¬ 
lattices of A sites and B sites, (a) A 2-ns-long spin-polarized current 
pulse (j = 3 x 10 13 A m -2 ) perpendicularly injected into an AFM 
nanodisk with a diameter of 80 nm in a circle region with a diam¬ 
eter of 40 nm followed by a 1-ns-long relaxation, (b) Close-up of 
the time evolution shown in (a), where the topological number of the 
two sublattices significantly changes from ~ 0 to d=l. 


spins continue to flip. As soon as we stop the current, an 
AFM skyrmion is nucleated to lower the DMI and AFM ex¬ 
change energies (see Supplementary Movie 1). It is relaxed 
to the optimized radius irrespective of the injected region, as 
shown in Fig. 2a (see Supplementary Movie 2). Once it is 
relaxed, it stays as it is for long, demonstrating its static sta¬ 
bility. We show the spin configuration of an AFM skyrmion 
obtained numerically in Fig. 2b. It is made of a toroidal DW 
with fixed radius and width determined by the material pa¬ 
rameters. There exists a threshold current density to create an 
AFM skyrmion, as shown in Fig. 2c. It is natural that the spins 
cannot be flipped if the injected current is not strong enough. 

The time-evolution of the topological charges of the AFM 
system is shown in Fig. 3. Note that there is a non-zero topo¬ 
logical number Q® in the AFM background state, which is 
created by the tilting magnetization at edges due to the DMI. 
It is Q A = —0.1827 for the A sites. The topological charge 
oscillates during the CPP injection. As soon as the CPP injec¬ 
tion is off, the topological charge develops suddenly to a fixed 
values. By subtracting Q° A from that of the AFM skyrmion in 
the A sites, we find Qa = 0.9865, which is almost 1. Simi¬ 
larly, we find Qb is almost —1. 

The AFM skyrmion can be created equally by a vertical 
current injection polarized along the -^-direction or the — z- 
direction (see Supplementary Movie 3). 

Phase diagrams. A skyrmion is topologically protected. 
Nevertheless, it may shrink or expand with the topological 
charge unchanged. We present a phase diagram in Fig. 4. It 
is convenient to understand it in terms of the DMI constant 
D. The DMI prevents a skyrmion from shrinking in antifer- 
romagnets as in the case of ferromagnets. (1) Near Dm 0, 
a skyrmion shrinks and disappears (blue region). (2) There 
are two cases when a skyrmion exists as a static stable object 
(yellow region): see also Fig. 2d. In one case (smaller D), 
its energy is more than that of the AFM ground state. It is 
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Figure 4. Phase diagrams of the AFM skyrmion as functions 
of the exchange stiffness J and the DMI constant D (a), and as 
functions of the PMA constant K and D (b). An AFM skyrmion 
cannot exist in the AFM ground (AFMG) state region (blue color). 
An AFM skyrmion can exist as a static stable object in the AFM 
skyrmion (AFMS) state region (yellow color). An AFM skyrmion is 
distorted to lower the DMI energy in the distorted AFM skyrmion (d- 
AFMS) state region (violet color). The ground state contains worm 
domains in the worm domain (WD) state region (green color). All 
spins are almost in plane because the exchange interaction or the 
anisotropy strength is too weak in the in-plane state region (light blue 
color), (c) Illustration of the states shown in the phase diagrams. The 
simulation is based on the nanodisk with a diameter of 80 nm. The 
J, D and K are in unit of 10 _21 J. The default values of J and K 
equal 15 x 10 -21 J and 0.8 x 10 -21 J, respectively. 


an energetically metastable state, but it is topologically sta¬ 
ble. In the other case (larger D), its energy is less than that 
of the AFM ground state. It would undergo condensation if it 
were not topologically protected. (3) When D becomes larger, 
a skyrmion is distorted to reduce the DMI energy, which we 
call a distorted AFM skyrmion (violet region). (4) When D 
becomes sufficiently large, a deformed skyrmion touches the 
edge, forming worm domains (green region). 

Creation of an AFM skyrmion from an AFM DW wall 
pair. A FM skyrmion can be created from a FM DW pair us¬ 
ing a junction geometry [11]. We show that a similar mecha¬ 
nism works in creating the AFM skyrmion as shown in Fig. 5. 
We first make an AFM DW pair through the CPP injection 
with p = —z. The AFM DW pair is shifted by applying a 
spin-polarized current through the STT on AFM DW [25] as 
shown in Fig. 5 (see the process from t = 10 ps to t = 20 
ps in Fig. 5). Here we consider the vertical injection of a 
spin current polarized along the —y- direction. The CPP in¬ 
jection moves the AFM DW in the rightward direction (+£). 
When the AFM DW arrives at the junction interface (t = 20 
ps), both the end spins of the DW are pinned at the junction, 
whereas the central part of the DW continues to move due to 
STT in the wide part of the nanotrack. Therefore, the structure 
is deformed into a curved shape and an AFM skyrmion tex¬ 
ture forms at t = 30 ps. This process is analogous to blowing 
air through soapy water using bubble pipes or plastic wands 
to create soap bubbles. The skyrmion will break away from 
the interface when the bulk of its structure continues to move 
rightward as shown at t = 40 ps. By continuously “blow- 
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Figure 5. Creation of an isolated AFM skyrmion in a 2D AFM 
nanotrack via AFM DW pair driven by a vertical spin-polarized 
current. A vertical current with a density of j = 4.5 x 10 11 A m -2 
in the wide part is applied perpendicular to the nanotrack from the 
bottom. An AFM skyrmion is created from an AFM DW pair driven 
by the current moving from the narrow part to the wide part of the 
nanotrack, D = 4 x 10 -21 J (see Supplementary Movie 4). The 
current density inside the wide part of the nanotrack is proportional 
to the current density inside the narrow part of the nanotrack with 
respect to the ratio of the narrow width (20 nm) to the wide width 
(100 nm). 


ing” AFM DWs into the junction, a train of AFM skyrmions 
is generated. 

Current-driven motion of an AFM skyrmion in a nan¬ 
otrack. We can move the AFM skyrmion by the CPP injection 
as in the case of the FM skyrmion. We show the relation be¬ 
tween the magnitude of the injected current and the velocity 
in Fig. 6 a, where the velocity is proportional to the injected 
current. 

We recall that the FM skyrmion is easily destroyed by 
touching the sample edges due to the SkHE. At the same time, 
the maximum velocity of the FM skyrmion in a FM nanotrack 
is typically much less than 10 2 ms -1 , limited by the confining 
force of ~ ( D 2 /J ) [26]. 

Conversely, there is no SkHE for the AFM skyrmion. 
Hence, it can move straightly in an AFM nanotrack without 
touching the edge. It is shown in Supplementary Movie 5, 
where a chain of encoded AFM skyrmions moves in a nan¬ 
otrack with a speed of ^ 1700 ms -1 driven by a vertical 
current without touching edges. 

In Fig. 6 b we compare the AFM and FM skyrmions. The 
velocity of AFM skyrmions can be very large compared to 
FM skyrmion, which is suitable for ultrafast information pro¬ 
cessing and communications. The steady motion of AFM 
skyrmions is demonstrated in Supplementary Movie 6 , where 
they move in a thin film without boundary effect driven by 
the vertical spin current. This highly contrasts with the case 
of FM skyrmions demonstrated in Supplementary Movie 7, 
where skyrmions do not move either parallel or perpendicular 
to the film edges. 

Discussion 

We have proposed magnetic skyrmions in the AFM system. 
The dynamics of AFM skyrmions is very different from those 
in the FM system, since they are topologically protected and 
are free from the SkHE. We have first checked that our simula¬ 
tion software reproduces a linear dispersion relation inherent 
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Figure 6. Motion of an isolated AFM (FM) skyrmion in a 
2D AFM (FM) nano track driven by a vertically injected spin- 
polarized current, (a) AFM skyrmion velocity and FM skyrmion 
velocity as functions of current density j with the CPP geometry. For 
the AFM skyrmion case, the current is applied along +z but polar¬ 
ized along —y,J = 15xlO _21 J. For the FM skyrmion case, the cur¬ 
rent is applied along +z but polarized along +y, J — — 15 x 10 -21 
J. For both cases, a = 0.3, D = 3.5 x 1(T 21 J ,K = 0.8 x 1(T 21 J. 
The open symbol denotes the destruction of the FM skyrmion due to 
the SkHE. (b) Top-views of vertical current-driven AFM skyrmion 
and FM skyrmion at selected current densities and times. 


to the two-sublattice structure of the AFM system, and then 
employ it to explore various properties of the AFM skyrmion. 
It is worth mentioning there are two recent preprints [27, 28] 
on AFM skyrmions, including the preliminary version of the 
present work [27]. Our work is focused on the injection 
and vertical spin current-driven dynamics of AFM skyrmions. 
Regarding the other work in Ref. 28, the thermal proper¬ 
ties as well as in-plane current-induced dynamics of an AFM 
skyrmion have been studied and a high-speed motion (v ~ 
10 3 m s -1 ) of an AFM skyrmion has also been shown in the 
absence of the SkHE, consistent with the present work. We 
believe that the AFM skyrmions will play a very significant 
role in the emerging field of AFM spintronics. 
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Methods 

Modeling and simulation. We perform the micromag- 
netic simulations using the Object Oriented MicroMagnetic 
Framework (OOMMF) together with the DMI extension mod¬ 
ule [14, 29-32]. The time-dependent magnetization dynamics 
is governed by the LLGS equation [33-37]. The OOMMF has 
been developed originally and used extensively for the sim¬ 
ulation of FM systems, and we have checked that one may 
use it to analyze the nanotexture in the AFM system as well. 
Indeed, we have successfully reproduced a linear dispersion 
relation inherent to the two-sublattice structure, as shown in 
Supplementary Figure 1. 

For micromagnetic simulations, we consider 0.4-nm-thick 
magnetic nanodisks and nanotracks on the substrate. With re¬ 
spect to the material parameters, we recall [38] that an anti- 
ferromagnet is a special case of a ferrimagnet for which both 
sublattices A and B have equal saturation magnetization. Both 
the DMI and the PM A arise from the spin orbit coupling, al¬ 
beit in different forms. We have checked that our results hold 


for a wide range of material parameters (cf. Fig. 4). Here, 
we use the parameters of the same order as those given in 
Ref. 39 for AFM materials. We thus adopt the magnetic pa¬ 
rameters from Refs. 6 and 14: the Gilbert damping coefficient 
a = 0.3, the gyromagnetic ratio 7 = —2.211 x 10 5 m A -1 
s -1 , the sublattice saturation magnetization Ms = 290 kA 
m -1 , the exchange constant J = 0 ^ 20 x 10 -21 J, the 
DMI constant D = 0^10xl0 - 21 J, and the PM A con¬ 
stant iC = 0^2xl0 _21 J unless otherwise specified. All 
samples are discretized into tetragonal cells of 1 nm x 1 nm x 
0.4 nm in the simulation, which ensures reasonable numerical 
accuracy as well as run time. The output time step is fixed at 
1 ps for the simulation of the dispersion relation, which is in¬ 
creased to 10 ps for the simulation of the skyrmion dynamics. 
The polarization rate of the spin-polarized current is defined 
as P = 0.4 in all simulations. The Zeeman field is set as zero 
because the AFM skyrmion, having no net magnetization, is 
insensitive to it (see Supplementary Figure 2). 
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